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Abstract 

The magnetic properties of a new family of molecular-based quasi- two dimension S =1/2 Heisen- 
berg antiferromagnets are reported. Three compounds, (Cu(pz)2 (0104)2, Cu(pz)2(BF4)2, and 
[Cu(pz)2(N03)](PFe)) contain similar planes of Cu 2+ ions linked into magnetically square lattices 
by bridging pyrazine molecules (pz = C4H4N2). The anions provide charge balance as well as 
isolation between the layers. Single crystal measurements of susceptibility and magnetization, as 
well as muon spin relaxation studies, reveal low ratios of Neel temperatures to exchange strengths 
(4.25/17.5 = 0.243, 3.80/15.3 = 0.248, and 3.05/10.8 = 0.282, respectively) while the ratio of the 
anisotropy fields Ha (kOe) to the saturation field Hsat (kOe) are small (2.6/490 = 5.3xl0~ 3 , 
2.4/430 = 5.5xl0 -3 , and 0.07/300 = 2.3xl0 -4 , respectively), demonstrating close approximations 
to a 2D Heisenberg model. The susceptibilities of CIO4 and BF4 show evidence of an exchange 
anisotropy crossover (Heisenberg to XY) at low temperatures; their ordering transitions are pri- 
marily driven by the XY behavior with the ultimate 3D transition appearing parasitically. The 
PFg compound remains Heisenb erg-like at all temperatures, with its transition to the Neel state 
due to the interlayer interactions. Effects of field-induced anisotropy have been observed. 

PACS numbers: 
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I. INTRODUCTION 

For more than seventy years |l| the study of low-dimensional magnetism has played 
an integral role in the understanding of phase transitions, critical behavior, and other 
aspects of quantum many-body physics. Highlights of this progression include the path- 
breaking neutron scattering studies of Jhe^ excitation spectrum of the one-dimensional 
(ID) S=l/2 Heisenberg antiferromagnet 
exchange-coupled layers of Cu(II) oxides 



the discovery of superconductivity in doped 
4j with the consequent flurry of theoretical |5|] and 
6|] research, as well as the discovery of macroscopic quantum tunneling in high- 



experimental t 
spin nanomagnets[3, Isj]. 

Experimental studies in quantum magnetism rely on the existence of simple and well- 
characterized model systems. For the case of the S — 1/2 (quantum) two-dimensional 
Heisenberg antiferromagnet on a square lattice (2D QHAF), the exchange must be as close 
to isotropic as possible and the layers must be well isolated. The appropriate Hamiltonian 
H, for a quasi-2D QHAF in an applied field H is 

H = JJ2 [S!S? + SfS] + (1 - A)S*Sj] 

nn 

+J' >^ Si - Sj> - ffAteH • Sj, (1) 

i,V j 

where the first summation is over nearest neighbors in the planes, the second summation 
links each spin to its counterparts in adjacent layers and the third includes all spins. Here J 
is the intralayer exchange parameter, J' is the coupling constant between spins in adjacent 
layers, and A is the exchange anisotropy parameter. (In this Hamiltonian, a positive J 
corresponds to antiferromagnetic exchange.) For an ideal 2D QHAF, J' = and A = 0. 
A full characterization of any physical realization of a model Hamiltonian must determine 
experimentally the values of J, J', and A. 

Little is known about the properties of the 2D QHAF in applied magnetic fields. The 
Zeeman term in the Hamiltonian (Eq. (TjQ)) has generally been neglected in studies of the 
copper oxides due to their large energy scale. Assuming the saturation field depends only 
on the exchange strength, the mean- field equation for H^at of an 5 = 1/2 system is given 
byfl 

Hsat = — — (2) 
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where z is the number of nearest neighbors and J is expressed in units of Kelvin. An 
exchange strength of 1500 K corresponds to an i^sAT ~ 4200 T and the usual experimental 
fields are minor perturbations to the Hamiltonian. 

Seven years ago some of the present authors published [lol. several studies of the 
molecular-based quasi-2D QHAF family which have exchange constants of 6.5 and 8.5 K, 
respectively, and respective saturation fields of 19 and 24 T. Although these materials 
did have relativel y l arge exchange ratios, J' /J ~ 0.24, they did show the characteris- 
tic susceptibilities [ 1 ol. 11] and heat capacities pjj of the 2D QHAF and revealed for the 
first time the characteristic upward curvature of the low-temperature magnetization of the 
2D QHAF. Preliminary reports of this discovery led to several theoretical studies of the 
magnetization [13I. \i\ and gave motivation for further investigations of the 2D QHAF in 



large applied fields. More recent papers 



15. 



161 ] predict an anomalous spin excitation spec- 



trum as the external field approaches H S at- In addition, it has been predicted [17|, ll8| that 
the application of an external field to a 2D QHAF would reduce quantum fluctuations along 
the field axis and gradually transform the spin anisotropy from Heisenberg toward XY, 
inducing a Berezinskii-Kosterlitz-Thouless (BKT) transition for perfectly isolated layers. 

In our previous work, we demonstrated the principle of using molecular-based magnetism 
to generate new families of 2D QHAF with moderate exchange strengths. In this paper, 
we report the magnetic properties of a new family of molecular-based quasi-2D QHAF 
(Cu(pz) 2 (C10 4 ) 2 , Cu(pz) 2 (BF 4 ) 2 , and [Cu(pz) 2 (N0 3 )](PF 6 )) in which Cu 2+ ions are linked 
into square magnetic lattices by pyrazine molecules (pz = C4H4N 2 ) with exchange strengths 
between 10 K and 18 K. The layers are well isolated by the counter-ions, as demonstrated 
by their T-^/J values of 0.24 to 0.30. The saturation fields range from 300 kOe to 490 kOe 
so available fields are powerful enough to test the recent theoretical predictions. (The 
three antiferromagnets under study will sometimes be identified as the C10 4 , BF 4 , and 
PF 6 compounds, respectively, in the interests of brevity.) 

The three compounds consist of Cu(II) ions bonded to four neutral bridging pyrazine 
molecules, creating positively charged 2D nets; the structures are charge-balanced by the 
counter-ions (ClOj, BFJ, NOj, and PFg ) which lie between the copper /pyrazine layers, 
with all anions except the PF 6 located in the axial sites of the copper atoms (Figs. [1] and 
[2]). A full description of the structures is found in Ref. [19. 



Each of three compounds forms magnetically square lattices, even though they have 
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three different low-temperature space groups. Cu(pz) 2 (C104)2 and Cu(pz) 2 (BF 4 ) 2 have C- 
centered monoclinic structures at low-temperatures (C2/c and C2/m, respectively), in which 
the copper sites are related by symmetry operations (c-glide in C2/c and C-centering in 
C2/m) that render every copper-copper bridge equivalent even though the crystallographic 
angles are 96.46° and 120.92°, respectively. The third compound, [Cu(pz) 2 (N0 3 )](PF 6 ), has 
a tetragonal space group 74 /mem in which the coppers sit on the four- fold rotation axis, 
generating square 2D copper/pyrazine nets. 

The copper/pyrazine layers are packed into three-dimensional lattices in two different 
ways, with important differences on the ultimate low-temperature 3D ordering transitions. 
For both Cu(pz) 2 (C10 4 ) 2 and Cu(pz) 2 (BF 4 ) 2 , the counter-ions are weakly coordinated to 
the axial sites of the copper atoms but do not bridge between the layers. To minimize steric 
hindrance, adjacent layers are offset by half a unit cell in along both axes within the layer so 
the counter- ions can interpenetrate, Fig. [T](b). Each copper atom is then equidistant from 
four coppers of the adjacent layer, an arrangement that leads to a cancellation of interlayer 
exchange interactions J'. 

The axial copper sites in the third compound, [Cu(pz) 2 (N03)](PF 6 ), are occupied by 
bridging nitrate groups (Fig. [2]) so there is a superexchange pathway between coppers of 
adjacent layers. For this reason, and the absence of any symmetry required cancellation 
of J', the interlayer coupling for [Cu(pz) 2 (N0 3 )](PF 6 ) should be the strongest of the three 
compounds. 



II. EXPERIMENTAL 

Crystals of Cu(pz) 2 X 2 were grown from aqueous solutions containing a 1:2 molar ratio of 
Cu(C104) 2 or Cu(BF 4 ) 2 and pyrazine. A drop of dilute HC10 4 (aq) or HBF 4 (aq) was added 
to the solutions to prevent precipitation of Cu(OH) 2 . Crystals of [Cu(pz) 2 (NOs)](PF 6 ) were 
grown from an aqueous solution of Cu(NOs) 2 with two equivalents of pyrazine and a 5-fold 
excess of KPF 6 after slow evaporation over several weeks. All crystals have the forms of 



dark blue tablets. Full details may be found elsewhere 

All DC magnetic susceptibility data were collected using Quantum Design MPMSR2 and 
MPMS-XL SQUID magnetometers. Susceptibility data on crystals of all three compounds 
were collected along all three orientation at various fields from 1.8 K to 60 K. High field mag- 
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FIG. 1: (a) Single layer structure of Cu(pz)2(ClC>4)2 at 293 K, viewed perpendicular to the ab 
plane, (b) Staggered layer structure of Cu(pz)2(C104)2- The CIO4 ions have been removed from 
the figure for clarity. 



netization data were collected on powder samples using a vibrating sample magnetometer 
(VSM) at the National High Field Magnet Laboratory in Tallahassee, Florida, in fields up to 
30 T at several temperatures. Because these fields were insufficient to saturate the magneti- 
zation data were collected in pulsed magnetic field experiments up to 60 T at the National 
High Field Magnet Laboratory at LANL. Zero field muon-spin relaxation (ZF /i + SR) mea- 



surements 



20 



2 If ] were made on powder samples of Cu(pz) 2 (BF 4 )2 and [Cu(pz) 2 (N03)](PF 6 ) 
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FIG. 2: Staggered layer structure of [Cu(pz) 2 (N0 3 )](PF 6 ) 



using the MuSR instrument at the ISIS facility, Rutherford Appleton Laboratory, UK and 
using the General Purpose Surface-Muon Instrument at the Swiss Muon Source, Paul Scher- 
rer Institute, CH. For these measurements samples were packed in Ag foil (thickness 25 /xm) 
and mounted on a Ag backing plate. 



III. RESULTS 



A. Powder susceptibility 

The molar magnetic susceptibilities (Xm) as a function of temperature for powder samples 
of Cu(pz) 2 (C10 4 ) 2 , Cu(pz) 2 (BF 4 ) 2 , and [Cu(pz) 2 (N0 3 )](PF 6 ) are shown in Fig. El These 
data were measured in a field of 1 kOe. For each compound, a broad rounded maximum is ob- 
served with the maximum value in Xm occurring near 15.9, 14.0, and 9.3 K for the CIO4, BF4, 
and [Cu(pz) 2 (N0 3 )](PF 6 ) compounds respectively. The data were compared to a numerical 
expression for the temperature dependent susceptibility of a 2D QHAF, as determined from 



high temperature series expansions [221] and quantum Monte Carlo simulations |23j] to obtain 



19 



exchange strengths and Curie constants as described in Ref 

The solid lines are the results of the best fits of the model to the data with fitting 
parameters (J, C) of (17.5(3) K, 0.426(6) cm 3 -K/mol), (15.3(3) K, 0.426(6) cm 3 -K/mol), 
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FIG. 3: Powder susceptibility measured in a 1 kOe field for Cu(pz)2 (0104)2 (o), Cu(pz)2(BF4)2 
(x) and [Cu(pz) 2 (N0 3 )](PF 6 ) (□). Solid lines represent fits to the 2D QHAF model. 



and (10.8(3) K, 0.439(6) cm 3 -K/mol) for the three compounds respectively. The quality 
of the fits is excellent except at the lowest temperatures where the data rise above the 
predicted susceptibilities. These discrepancies are discussed in detail below. The value of 
the Curie constant for the perchlorate compound corresponds to an average ^-factor of 2.13, 
in excellent agreement with the room temperature value found from the EPR measurements. 
Fuller discussion of the powder susceptibilities and the magnetostructural correlations in the 



copper pyrazine family are found elsewhere 



19J. 



B. High field magnetization 

Fig. H](a) shows the relative magnetization M(H)/Msat as a function of field for poly- 
crystalline samples of all three copper pyrazine compounds. Data were collected for each 
compound at several temperatures between T=0.50 K and 4.00 K with representative points 
of the the T=0.5 K data shown in Figs. H^a) andH^b). For every compound the behavior 
is similar: a small region in which the magnetization is linear in field, followed by a gradual 
upward curvature that continues to nearly the saturation value. Rounding due to thermal 
excitations and the variation in g-f actor values is negligible at 0.5 K for all three samples 
as seen in the PF6 data set. The observed rounding for the BF4 and CIO4 data near their 
saturation fields is due to their reaching Msat near the peak field when dH/dt is small, a 
feature in pulsed field magnetization. Values for the saturation fields were estimated to be 



S 



310(15) kOe, 430(20) kOe, and 520(20) kOe for the PF 6 , BF 4 , and C10 4 compounds respec- 
tively, taking into consideration the rounding of the BF4, and CIO4 data sets curves near 
saturation. The zero-temperature saturation field of a 2D QHAF with exchange strength 




FIG. 4: (a) M/M S at vs. fi H at 0.50 K for Cu(pz) 2 (C10 4 ) 2 (O), Cu(pz) 2 (BF 4 )2 (□) and 
[Cu(pz) 2 (N0 3 )](PF 6 ) (A). 

(b) M/Msat vs. H/H SAT at 0.50 K for Cu(pz) 2 (C10 4 ) 2 (O), Cu(pz) 2 (BF 4 ) 2 (□), and 



2J], the solid 



[Cu(pz) 2 (N0 3 )](PF 6 ) (A). The short dash line is the T = ID QHAF model 
black line is a Monte Carlo Calculation of a 2D QHAF at T/J = 0.05(Ref l25l). and the solid gray 
curve is a T = spinwave calculation [141 ] . 



J and four nearest neighbors can be calculated using Eq. (j2J). Based on the values of the 
exchange strengths determined by the powder susceptibility measurements and using an 
average g- value of 2.13 for each compound, the predicted saturation fields for the three 
compounds are 302(8), 428(8), and 489(8) kOe, respectively, with the uncertainty of the 
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calculated saturation fields due to the uncertainty of the experimental exchange constants. 
Within these uncertainties, the predicted T = saturation fields are equal to the experi- 
mental values obtained at 0.5 K for the PF 6 and BF 4 compounds, but about 6 percent lower 
than the experimental value of the CIO4 data set. Given the lower quality of the high-field 
data for the CIO4 compound, the value of Hsat = 490 kOe obtained from Eq. §Z§ will be 
used henceforth while values of 300(10) kOe and 430(10) kOe will be used for the PF 6 and 
BF4 compounds. In order to compare the three data sets to each other, the relative mag- 
netizations were plotted as a function of the relative fields if/if sat in Fig. IU(b). As seen in 
the figure, near universal behavior occurs. 

Included in Fig. IH(b) are results from numerical calculations of M(H)(T = 0) based 
on spin wave calculations for a 2D QHAfQ solid gray line, and a Monte Carlo 



simulation 



25| of a 2D QHAF finite size lattice at T/J = 0.05 (solid black line). The 



two lines are indistinguishable except just below the saturation field. Also included is a 
short dash line representing the T = magnetization of a one-dimensional QHAF. See the 
discussion in section HVl for details. 



C. Single Crystal low field Magnetization and Susceptibility 

Cu(pz) 2 (C104)2- Crystals of Cu(pz) 2 (C104)2 grow as flat plates with the cop- 
per / pyrazine layers lying in the fee-planes, parallel to the dominant face Ji| . The axis normal 
to the layers is a* which will be defined as the ^-direction. The two magnetically equiva- 
lent axes in the layer are defined as x. The low field magnetization of single crystal of the 
perchlorate salt was studied at low temperatures along three orthogonal directions, two in 
the plane of the layers and one normal to the layers. The response of the two directions in 
the layers were identical, even through the crystal is monoclinic (C2/c at low temperatures) 
and not tetragonal. The behavior normal to the layers is different. 

Fig. (jSj) shows the molar magnetization of Cu(pz) 2 (C10 4 )2 at 1.8 K between zero and 
10 kOe with the fields aligned in the plane (x) and normal to the plane (o). When the field 
is perpendicular to the layers, the magnetization increases linearly throughout the displayed 
range. Only for fields exceeding several tesla does the slope gradually begin to increase and 
display the high field behavior seen in Fig. H](a). In contrast, for fields within the layers, 
the magnetization shows two linear regions with the smaller slope occurring between and 
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4000 6000 
H(Oe) 



FIG. 5: Low field magnetization data of Cu(pz)2(C104)2 for H || x (x) and H || z (o) at T = 1.8 K. 



2.6(5) kOe. No further breaks of slope occur at higher fields, only a gradual increase before 
saturation. 

The two susceptibilities, Xx and Xz, as determined in a field of 1 kOe, are shown at low 
temperatures in Fig. [61 Two types of anisotropy are observed, with the first being the simple 
g-factor anisotropy g z = 2.27 and g x = 2.07 as determined by EPR. The corrected suscep- 
tibilities, Xi/Ci, are shown in the lower panel of the figure where the second, temperature 
dependent, anisotropy is evident. As the temperature drops below that of the susceptibility 
maximum, Xx decreases steadily, though not to zero, while a minimum is found in Xz at 
4.5 K; for temperatures below the minimum, the susceptibility increases linearly upward. 
Zero-field AC-susceptibility measurements confirm the 1 kOe DC susceptibility presented in 
Fig [6j As will be discussed below, the minimum is neither the signature of long-range order 
nor a field- induced effect, but is due to an intrinsic spin-anisotropy crossover. 

Cu(pz) 2 (BF 4 ) 2 . The behavior of Cu(pz) 2 (BF 4 ) 2 at low fields is very similar to 
Cu(pz) 2 (C104) 2 , with the magnetization curve showing a slope change at around 2.5(6) kOe 
for fields within the layer while no such breaks were observed for a field perpendicular to 
the layer. The M/H vs. T data also show a minimum around 4.05 K only when the field is 
parallel to the layer. 

[Cu(pz) 2 (N03)](PF@). This tetragonal compound also grows as flat plates. The normal 
to the plate is the c-axis and will be referred to as the z-direction. The equivalent axes in 
the plane are the tetragonal a-axes and point in the x-direction. The magnetization as a 
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T(K) 



FIG. 6: (a) Susceptibility of Cu(pz) 2 (C10 4 ) 2 at 1 kOe for both H \\ x (x) and H \\ z (o). Below 
5 K Xz rises above and Xx falls below the solid lines representing the susceptibility for an ideal 2D 
QHAF. The inset shows the details of Xz between 2 K and 6 K. (b) Susceptibility data normalized 
by Curie constant for Cu(pz)2(C104)2, same notation as in (a). 



function of field of the PF 6 compound at 1.8 K is shown in Fig. [71 A field-induced transition 
at around 70 Oe is observed for a field along x axis takes place while no such anomaly 
occurs for fields in the z direction. This behavior is similar to that seen for the CIO4 and 
BF4 compounds but at a substantially smaller field. For fields greater than 100 Oe the 
magnetization along each axis is linear up to 5 kOe, then gradually curves upward. M z is 
greater than M x for the same field and temperature, consistent with the greater value of g z . 

The zero-field ac susceptibilities for [Cu(pz) 2 (N0 3 )](PF 6 ) (filled symbols in Fig. [81(a) have 
similar temperature dependences as those of both the the zero-field ac susceptibilities and 
DC susceptibilities of Cu(pz) 2 (004)2 (FigEJ; there is a minimum for \ z but no minimum in 
Xx- However, for [Cu(pz) 2 (N03)](PF 6 ) in DC fields of 100 Oe or greater, a minimum appears 
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FIG. 7: Single crystal magnetization data (T = 1.8 K) of [Cu(pz) 2 (N0 3 )](PF 6 ) crystal for H \\ x 
(x) and H \\ z (o). 



for Xx as well. Data measured in a field of 1 kOe are shown in Fig. [HJ The mimimum for Xz 
is intrinsic and due to an internal anisotropy, as is the case for Cu(pz)2(C10 4 ) 2 ; the addition 
of the DC-field modifies the minimum but does not create it. In contrast, the minimum for 
Xx in the 1 kOe field is not intrinsic, but is field-induced. An analysis of the two types of 
anisotropies is found in the Discussion. 



D. Determination of Tn 



The unambiguous identification of long range magnetic order in low-dimensional systems 
is often made difficult by the existence of quantum fluctuations, which act to depress the 
ordering temperature and reduce the size of the magnetic moments. This often hinders 
experimental methods such as magnetic susceptibility or magnetic neutron diffraction. In 
addition, the build up of spin correlations above the ordering temperature in these systems 
reduces the available entropy and hence the response of heat capacity 2lj . In contrast, 
muon-spin relaxation measurements have been shown to detect magnetic order in cases 
where transitions are very difficult to observe with more conventional techniques 26] . 

Cu(pz) 2 (C104)2 The zero-field Neel temperature of Cu(pz) 2 (C104)2 has recently been 
determined to be 4.21(1) K using /i + SR 21]. Using single crystal magnetometry, we have 
confirmed this value for the CIO4 compounds and have also precisely ascertained the ordering 
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FIG. 8: (a) M/H ratios of [Cu(pz) 2 (N0 3 )](PF 6 ) for both H \\ x (A) and H \\ z (o). The inset 
shows details of M/H for both directions between 2 K and 6 K. Open symbols represent DC data at 
1 kOe and filled symbols represent AC data measured at zero-field, (b) DC M/H data normalized 
by Curie constant for [Cu(pz)2(N03)](PFe), same notation as in (a). 



temperatures for the BF 4 and PF 6 compounds. 

In Fig. [HI starting at 1.8 K, for 1 kOe, the susceptibility along z-direction drops as 
the temperature rises until, at a temperature just below that of the minimum, the slope 
of the curve drops sharply at 4.25 K, before undergoing a change of curvature with the 
slope growing more positive (see inset, Fig. Eta)). The sudden change of curvature occurs 
at the same temperature at which the muon experiment 2l| found the internal fields of 
Cu(pz)2 (0104)2 to vanish, within experimental error. 

Cu(pz) 2 (BF 4 ) 2 . ZF /i + SR spectra measured on Cu(pz) 2 (BF 4 ) 2 at several temperatures 
are shown in Fig. (9^a). Below Tn (Fig. Ws)) we observe oscillations in the time dependence 
of the muon polarization (the "asymmetry" A(t) ) which are characteristic of a quasi-static 
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FIG. 9: (a) Temperature evolution of ZF /i + SR spectra measured on Cu(pz)2BF4. (b) Above Tn 
low frequency oscillations are observed due to the dipole-dipole coupling of F-^t + -F states. Inset: 
The energy level structure allows three transitions, leading to three observed frequencies, (c) Below 
Tn higher frequency oscillations are observed due to quasistatic magnetic fields at the muon sites. 
Inset: The evolution of the muon precession frequency Va with temperature. 



local magnetic field at the muon stopping site. This local field causes a coherent precession of 
the spins of those muons for which a component of their spin polarization lies perpendicular 
to this local field (expected to be 2/3 of the total spin polarization for a powder sample). The 
frequency of the oscillations is given by z/j = 7 M |2?j|/27r, where 7 M is the muon gyromagnetic 
ratio (= 27r x 135.5 MHz T -1 ) and Bi is the average magnitude of the local magnetic field 
at the ith muon site. The precession frequency v is proportional to the order parameter of 
the system. Any fluctuation in magnitude of these local fields will result in a relaxation of 
the oscillating signal, described by relaxation rates \. The presence of these oscillations at 
low temperatures in Cu(pz)2(BF4)2 suggests very strongly that this material is magnetically 
ordered below T N = 3.7(1) K. Oscillations are observed at a single frequency, suggesting a 
single muon site in this material (in contrast to the case of the CIO4 material, where three 
frequencies were observed 21J] ) . We note, however, that the oscillations are quite heavily 
damped (with a typical relaxation rate of A ~ 1.7 MHz) compared with those measured in 
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the CIO4 case (where A ~ 0.4 MHz) which may imply the presence of several muon sites 
experiencing slightly different local magnetic fields. 

Above T N the character of the measured spectra changes considerably (FigJS^a) and (b)) 
and we observe lower frequency oscillations characteristic of the dipole-dipole interaction of 
the muon and the 19 F nuclei comprising the BF4 counter-ions. This behavior 27j has been 
observed previously in systems of this kind 28j]. Fits to the data show that the muon stops 
between two fluorine atoms and the resulting F-/i + -F spin system consists of four distinct 
energy levels with three allowed transitions between them (inset, Fig.[9](b)) giving rise to the 
distinctive three-frequency oscillations observed. Fits of the data suggest a fi + -F separation 
of 0.12 nm and an F-/i + -F bond angle of ~ 140°. 

Variable temperature measurements of the magnetization of single crystals of 
Cu(pz)2(BF4)2 (not shown) have revealed behavior similar to that seen in the CIO4 com- 
pound shown in Fig. [6j although at slightly lower temperatures. A slope anomaly indicative 
of the ordering temperature and the temperature of Xmin of Cu(pz) 2 (BF 4 ) 2 have been found 
to be 3.80(5) K and 4.05(5) K, respectively. 

[Cu(pz) 2 (N03)](PF 6 ). Variable temperature studies of single crystals of the PF 6 
compound (inset of Fig. [8]) show that the ordering temperature and the temperature of 
Xmin are very close to each other, both occurring at 3.05 K. The ordering temperature of 
[Cu(pz) 2 (N0 3 )](PF 6 ), as determined by /i + SR experiments, has been previously reported}^ 
to be near 2.0(2) K. The significant disagreement between the two values for the critical tem- 
perature of [Cu(pz) 2 (N03)](PFg) is in sharp contrast to the excellent agreement found from 
the two techniques for Cu(pz) 2 (C104) 2 and the good agreement found for Cu(pz) 2 (BF 4 ) 2 . It 
is worth noting, however, that the oscillations observed for this material due the the magnetic 
ordering were very heavily damped (with A ~ 4 MHz) , making the unambiguous determina- 
tion of an ordering temperature very difficult in this case. The previous estimate obtained 
from the muon measurements, corresponding to the disappearance of the magnetic oscilla- 
tions (due to the relaxation rate increases with increasing temperature), should therefore be 
taken as a lower bound on Tn. The heavy damping observed in this material is suggestive of 
a significantly larger width of the local magnetic field distribution in [Cu(pz) 2 (NC>3)](PF 6 ) 
compared to the CIO4 and PF 6 materials. Our yU + SR measurements do, however, unambigu- 
ously confirm the presence of long range magnetic order in this material at low temperature. 
For the purposes of this paper, we will consider the Neel temperatures to correspond to the 
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temperatures of the slope anomalies in \ z . 



IV. DISCUSSION 

A. High Field Magnetism 

The T = 0.5 K magnetizations of CIO4, BF 4 , and PF 6 materials are qualitatively similar; 
They begin with a small initial slope that gradually increases with field until the maximum 
slope is reached just before saturation; the respective saturation fields as based on the 
exchange constants of the three compounds are calculated to be 490, 430, and 300 kOe 
(Fig. H^a) ). When appropriately normalized as H/ H$at, the three compounds show universal 



behavior (Fig. H(b)). Similar behavior was observed in a previous study lOj, llll] of several 
less well-isolated copper bromide 2D QHAF. 

The T = magnetization for the 2D QHAF has been calculated using a spin-wave 
expansion [141] and is displayed as the solid gray curve in Fig.0](b). It is qualitatively similar to 
the normalized data but has a somewhat smaller slope for the first half of the magnetization 
curve and rises to an infinite slope at the saturation field. In addition, the magnetization 
curve at a relative temperature k^T/J = 0.05 has been calculated using a quantum Monte 
Carlo simulation [23J] and plotted as the solid black line in Fig.|U(b). (This simulated result is 
identical to the T = calculation except for values of H/ Hsat > 0.9 where is slightly reduced 
due to the thermal excitation of spin waves.) The discrepancy between the experimental and 
theoretical slopes is attributed to the fact that all three compounds are in the 3D ordered 
state at T = 0.5 K while the model describes an ideal 2D QHAF. As seen in Fig. [3], the 
experimental low-temperatures susceptibilities are all higher than those predicted by the 2D 
QHAF model. The magnetization data of these 2D compounds display much less curvature 
that appears in the ID QHAF, represented by the short dash line in Fig. IU(b); the difference 
arises from the greater influence of quantum fluctuations in the ID system which reduces 
the effective moments. 

High-field magnetization of a similar series of copper /pyrazine 2D antiferromagnets 



(]Cu(pz) 2 (HF 2 )](X), where X =C10 4 , BF 4 , PF 6 , SbF 6 ) have recently [29fl been reported 
on both single crystal and polycrystalline samples. These results for the polycrystalline 
samples are similar to those displayed in Fig. H] but the magnetization curves for single 
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TABLE I: Parameters for the layered compounds Cu(pz)2(C104)2, Cu(pz)2(BF4)2, 
[Cu(pz) 2 (N0 3 )](PF 6 ) and Sr 2 Cu0 2 Cl2. 





J 


7n 




J'/ J a 


H A 


-HsAT 


-Ha/-£^sat 


Aco 






(K) 


(K) 






(kOe) 


(kOe) 








Cu(pz) 2 (C10 4 )2 


17.5 


4.25 


0.243 


8.8 xl0~ 4 


2.6 


490 


5.3xl0~ 3 


4.6x10" 


-3 


Cu(pz) 2 (BF 4 )2 


15.3 


3.8 


0.248 


l.lxlO" 3 


2.5 


430 


5.8xl0~ 3 


6.2x10" 


-3 


[Cu(pz) 2 (N0 3 )](PF 6 ) 


10.8 


3.05 


0.282 


3.3xl0~ 3 


0.07 


300 


2.3xl0~ 4 


1.2x10" 


-2 


Sr 2 Cu0 2 Cl2 


1450 


251 


0.173 


1.9xl0" 5 


7 


40000 


1.8xl0" 4 


3.6x10" 


-4 



"Based on Eq. (HJ which assumes A = 



crystals are much less rounded near the saturation field, a consequence of the presence of a 
single g-factor. For typical copper compounds, the ^-factors range from 2.07 to 2.27, leading 
to a spread of critical fields of roughly ten percent in studies of randomly oriented samples. 



B. Interlayer Exchange and Intrinsic Anisotropics 



A figure of merit of characterizing low- dimensional antiferromagnets is the critical ratio 
of the Neel temperature to the dominant exchange strength, T N /J. For Cu(pz) 2 (C104) 2 , 
Cu(pz) 2 (BF 4 ) 2 , and [Cu(pz) 2 (N0 3 )](PF 6 ), the critical ratios are 0.24, 0.25 and 0.28, respec- 
tively (Tabled]). The ideal 2D QHAF has a critical ratio of zero 30j so the ordering present in 
these compounds must arise from lattice anisotropy (3D interaction J') , exchange anisotropy 
(A 7^ in the exchange Hamiltonian in Eq. (OQ)) or both. As will be shown below, it is not 
possible to determine an accurate value of J' in the presence of a nonzero A. 

An alternative measure of the degree of isolation for quasi-2D QHAF is given by the 
magnetic correlation length £ at the ordering temperature £(Tn)). The correlation length 
is known to diverge exponentially at low temperature with only a weak temperature de- 



pendence in the prefactor 



3l| 



lattice constant a is given by 32 ] 



'he full expression for the correlation length in units of the 



£ e c/a ( 2np 

~ = qo ex P 

a 8 2rcp s \ 1 



1-0.5- 



T 



+ 0( 



T 



(3) 



where c = 1.657 J a and p s 



2ixp s 2-Kps 
0.1830 J are the renormalized spin wave velocity and spin 
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stiffness constants, respectively. For Cu(pz) 2 (0104)2 with a critical ratio of 0.24, the cor- 
relation length at the ordering temperature is predicted to be £(Tn)/o = 50, while the 
values for Cu(pz) 2 (BF 4 ) 2 , and [Cu(pz) 2 (N0 3 )](PF 6 ) are calculated to be 45 and 25, respec- 
tively. For comparison, the correlation lengths £/a at Tn for Sr 2 Cu0 2 Cl 2 and deuterated 
copper formate tetrahydrate have been determined by neutron scattering experiments to be 
approximately 220 (Ref. 33) and 55 (Ref. 3^), respectively. 

The ideal 2D QHAF only orders at T = 0, but all quasi-2D QHAF contain finite 
interlayer couplings J' that induce long-range order (LRO) at a finite temperature T N . 
Recently a method for estimating the interlayer coupling constant J' in 3D arrays of isotropic 
2D QHAF has been developed based on a modified random phase approximation, modeled 



with classical and quantum Monte Carlo simulations 
formula relating J' and Tn, 



351 ] . The approach leads to an empirical 



J ' «pfl-^l, (4) 



J \ 2~N / 

where p s is the spin stiffness (p s = 0.183 J for the 2D QHAF [36]) and b = 2.43 for S = 1/2. 
This result shows that T^/J decreases only logarithmically with J'/J ratio; even very well 
isolated 2D layers will have critical ratios far from zero. For J 1 j J ratios of 10 _1 , 10~ 2 , 10~ 3 , 
the corresponding values of T^/J are 0.491, 0.326, and 0.244. The 2D QHAF known to 
have the greatest degree of isolation is Sr 2 Cu0 2 Cl 2 with a Tn/J = 256.5 K/1450 K = 0.18 
(Ref. m), corresponding to J'/J=3xl()- 5 . For La 2 Cu0 4 , T N /J = 311 K/1500 K = 0.207, 
corresponding to J'/J= 1.7xl0 -4 . Using Eq. and the Tn/J ratios for the CIO4, BF 4 
and PF 6 compounds, their J'/J ratios are found to be 8.8xl0~ 4 , l.lxlO -3 , and 3.3 xlO -3 
respectively. 

The values of J' obtained from Eq. (J4]) are not necessarily correct for a given compound, 
since they are calculated under the assumption that the exchange interactions are strictly 
Heisenberg, A = 0. Purely Heisenberg magnetic behavior can be found in systems of organic 
radicals or transition metal ions with spin-only moments, due either to half filled shells or 
complete quenching of orbital angular momentum. The 3d 5 Mn(II) has a half-filled shell 
and a spin-only moment, as ascertained by its ^-factor of 2.00. However, the orbital angular 
momentum of the 3d 9 Cu(II) ion in a non-cubic site is not completely quenched; enough 
orbital angular momentum remains to create ^-factors ranging typically from 2.05 to 2.25 
for different orientations of the external field relative to the coordination geometry. These 
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remaining internal fields at the copper site typically lead to anisotropy parameters in the 
range -0.02 < A < 0.02. 

Evidence for exchange anisotropy is found in the low-temperature, low-field magnetization 
curves (Figs. [5] and [7]). For all three compounds, changes of slope in the magnetization curve 
are found when the field is applied within the copper /pyrazine planes; no magnetization 
anomalies occur for fields normal to the planes. These anisotropy fields, labeled as Ha, are 
2.6 kOe, 2.5 kOe, and 70 Oe for the CIO4, BF 4 , and PF 6 compounds, respectively, Table fl] 
Normalized by their respective saturation fields -£/sat of ~ 490 kOe, 430 kOe, and 300 kOe, 
the anisotropy ratios H a /H S at, are 5.3 x 10" 3 , 5.8 x 10" 3 , and 2.3 x 10~ 4 (Tabled]). The 
ratios for CIO4 and BF 4 are close to a percent but the value for the PF 6 is surprisingly more 
than a factor of thirty smaller. These anomalies point to an XY anisotropy as they occur 
for fields applied within copper pyrazine planes, normal to the anisotropy axis; no anomalies 
are found for fields along the z-ax.es. 

Additional evidence for XY anisotropy is found in the single crystal susceptibilities. 
Were Ising anisotropy present, the easy axis susceptibility would descend to zero as the 
temperature drops below the Neel temperature. In the presence of an XY anisotropy, all 
susceptibilities remain finite in the zero-temperature limit due to the continuous rotational 
symmetry of the ground state in the xy plane; it is this behavior that is observed for all 
measured susceptibilities. Yet it is the minima in the susceptibilities \ z for each compound 
that provide the clearest evidence of XY-anisotropy. 

The minimum in the out-of-plane component of the uniform susceptibility of a 2D QHAF 



has previously been recognized as a signature of XY anisotropy [13, llB] . As discussed below, 
the temperature of the minimum in Xz, denoted as Tco, marks the crossover from Heisenberg 
to XY behavior. Based on the ratio of Tco to the intralayer exchange constant J, Cuccoli 
and coworkers have obtained an empirical formula [18J] which can be used obtain an estimate 
for the exchange anisotropy parameter Aco- 

4vr ■ 0.214J 



-co 



(5) 



ln(160/A co ) 

Implementation of that formula for the CIO4, BF 4 , and PF 6 compounds leads to values 
A co which are 4.6xl0" 3 , 6.2 x 10 -3 , and 1.2 x 10~ 2 (Tabled]). We note that these values 
for the anisotropy parameters for the CIO4 and BF 4 compounds are very close to the values 
of the respective field ratios Ha/ Hsat and conclude that the minimum in Xz does provide a 
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quantitative measure of the degree of XF-anisotropy in quasi-2D QHAF. 

For the PF 6 compound the values of H a /H SA t (2.3x10~ 4 ) and A co (1.2xl(T 2 ) differ by 
a factor of 50. The compound with the smallest anisotropy has the highest relative crossover 
and ordering temperatures. It is likely that this discrepancy is due to a high value of J' / J 
for the PF 6 compound. Eq. §5§ was derived on the assumption that J'— so the existence 
of a minimum in Xz would only be due to an exchange anisotropy. For the PF$ compound 
the 3D ordering temperature dominates and prevents the material from ever reaching a low 
enough temperature for the exchange crossover to appear. 

There is a simple qualitative explanation for the appearance of the minimum in Xz in a 
magnetic system with XY anisotropy. At high temperatures the system is isotropic with 
the orientation of antiferromagnetically coupled pairs fluctuating in all three directions. Xx 
and Xz are equal (within (/-factor anisotropy) and they decrease equally with decreasing 
temperature (T < J) as there is less thermal energy to overcome the antiferromagnetic 
coupling. As the temperature is lowered further, the XY anisotropy becomes increasingly 
relevant and an larger fraction of the spins anti-align in the plane. For a field in the z- 
direction, this process increases the fraction of antiferromagnetically coupled spins which 
will cant in the direction of the field; consequently, Xz begins to increase as the temperature 
cools and its minimum appears. In contrast, for a field in the plane, the number of responding 
spin pairs canting in the direction of the field decreases (there are fewer moments along +/- 
z) and Xx falls below the value of the isotropic 2D QHAF (See Fig. [6(a)). The minimum 
in Xz marks the temperature at which the out-of-plane component of the antiferromagnetic 
coupling becomes irrelevant and the system crosses over from Heisenberg to XY anisotropy. 

In light of the presence of the XY anisotropy, the physical significance of the anisotropy 
fields becomes apparent. At the low temperature of 1.8 K, the antiferromagnetically-coupled 
moments fluctuate within the copper /pyrazine planes. Application of a field within the plane 
will cant the moments oriented primarily normal to the field, but there will be little response 
of the moments along the axis of the field. As the field is increased towards if a the Zeeman 
energy overcomes the anisotropy energy and the moments are now free to rotate out of the 
easy plane. This dramatically increases the number of moments capable of canting in the 
direction of the field; the slope of the M vs. H curve should increase substantially for fields 
greater than H^. Experimentally the ratios of the M vs. H slopes above the critical field 
to the slopes below have been found to be 1.94, 1.68, and 1.37 for the CIO4, BF4 and PF@ 
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compounds, respectively. The small value for the PFg compound is consistent with its much 
smaller anisotropy parameter. 

A Xz minimum has previously been reported 37| for Sr 2 Cu0 2 Cl 2 , the 2D QHAF with 
the lowest critical ratio, (Table [I]). Magnetic susceptibility and 35 C1 nuclear spin-lattice 



relaxation rate studies 
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39] were interpreted as demonstrating the existence of a small 
amount of XY anisotropy. According to Eq. (jSJ), the ratio of crossover temperature to 
exchange strength for this compound (300 K/1450 K = 0.21) corresponds to A<x> = 3.6 x 
10~ 4 . Measurements of the out-of-plane spin wave gap by neutron scattering 33] directly 
established the value of the anisotropy parameter to be A cxp = 1.4 x 10~ 4 , in very good 
agreement with the formula. Similar experiments are required to directly determine the 
anisotropy parameters of the copper/pyrazine antiferromagnets. 



C. Long range order 



In the presence of exchange anisotropy, two-dimensional magnets order at finite temper- 
atures even in the absence of 3D interactions. If A < 0, the S x and S y components are 



more heavily weighted than S z ; in the limit A = 1, the spin degrees o: 
totally reduced from three to two and spontaneous order occurs 
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freedom have been 



41] at the Berezinskii- 



Kosterlitz-Thouless (BKT) transition temperature, Tbkt = 0.90JS 2 for classical spins and 
0.353 J for 5 = 1/2 (Ref |42j). This unique order is characterized by a diverging suscepti- 
bility but no spontaneous magnetization. No magnetic system is known to undergo a such 
a transition but BKT behavior has been observed in superfluid or superconducting films, 
as well as in arrays of Josephson junctions [43]. Similarly, for Ising-like anisotropy (A < 0), 
the increased weighting of the axial components leads to spin-spin correlations diverging at 
finite temperatures. In the presence of complete Ising anisotropy (A — > oo), Tq = 2. 269 J 
(Ref |44|). Given the importance of exchange anisotropy in determining the critical ratio, no 
estimate of J' is possible without first determining the sign and strength of A. We have done 
so for Cu(pz)2(C10 4 ) 2 , Cu(pz) 2 (BF 4 ) 2 , and [Cu(pz) 2 (N0 3 )](PF 6 ) in the previous section. 

Can the presence of a small XY anisotropy (~10 -2 ) induce long-range order? While pure 
XY ansiotropy (A = 1 in Eq. (TjQ)) leads to the BKT transition, quantum fluctuations have 
a larger influence near the Heisenberg limit of A = 0. It was not known that small XY 
anisotropies were capable of stablizing order until the QMC studies of Ding 45]. Simulations 
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on lattices up to 96 x 96 sites demonstrated that BKT order is surprisingly stable in the 
presence of small anisotropies. For anisotropy values of 0.02, 0.1, and 0.5, the critical 
temperatures were found to be 0.250 J, 0.285J, and 0.325J, respectively. Scaling arguments 
show that the naive guess that the critical temperatures should scale with the anisotropy 
parameter is incorrect; instead, Tbkt ~ l/ln(C/A). More recent simulations on larger 



lattices [17j confirmed the conclusions of Ding and extended the study to an anisotropy as 
small as 10~ 3 , finding that Tbkt = 0.175 J and 0.229 J for anisotropies of 10~ 3 and 0.02. The 
dependence of their critical temperatures upon the anisotropy parameter can be expressed 
by the second empirical equation 

2 22 

TbKt/ J = ln(330/A) (6) 

This equation is plotted as the solid line in Fig. [TUJ. It is based on the assumption 
that J' = 0; the presence of a finite J' would raise the curve. By putting in the values 
of the exchange anisotropy for the three copper/pyrazine compounds as determined by the 
Hj±/H S at, Table [H Eq. (jSJ) can be used to predict the critical ratios for the three compounds, 
in the absence of the J' interaction. For the CIO4, BF 4 , and PF 6 compounds, the critical 
ratios would be 0.20, 0.22, and 0.16, respectively. The corresponding value for Sr 2 Cu0 2 Cl 2 is 
0.15. The experimental critical ratios of the same four compounds are also plotted on Fig.fTOl 
In this plot, if the transition is primarily driven by the XY anisotropy, one anticipates the 
experimental data will be found somewhat above the curve, the exact amount determined by 
J'. At the opposite limit, with the transition primarily driven by the interlayer interactions, 
the experimental value will appear substantially above the curve. 

Fig. [T0l shows three of the compounds follow a similar pattern while the fourth is a definite 
outlier. The values for C10 4 , BF 4 , and Sr 2 Cu0 2 Cl2 follow the general shape of the curve 
but lie above it by 0.04 (0.24 - 0.20) for C10 4 , 0.02 units for BF 4 and 0.03 for Sr 2 Cu0 2 Cl 2 . 
This behavior is consistent with their description as quasi-2D systems in which the stronger 
perturbation to the basic 2D QHAF Hamiltonian is the moderate (~ 0.5% for CIO4 and 
BF 4 ) XY anisotropy. In the absence of interlayer interactions, the critical ratios for these 
three compounds would have been been located on the solid curve but the presence of J' 
enhanced the transition temperature approximately 20% to their measured values. Since 
their ordering is primarily due to the anisotropy, the values for J' / J in Table |U(8 x 10~ 4 
and 1 x 10~ 3 , respectively) derived from Eq. (T4J) (which assumed A = 0) are clearly too 
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FIG. 10: The critical ratio Tn/J is plotted against the exchange anisotropy parameter A = 0. The 



17 ]. The four experimental values are 



solid line is the prediction of Eq. ([6]) which assumes J' = 
based on the information given in Table [J with A = Ha/Hsat- The error bars of the data are 
smaller than the sizes of the symbol. 



high. The appropriate values of J'/ J are those which will raise the critical ratio of 0.20 for 
a purely 2D anisotropic magnet to 0.24 for a quasi-2D anisotropic magnet. In the absence 
of guidance by simulations, it is not useful to speculate about the appropriate values. We 
do note that the closer proximity of the BF 4 value to the theoretical curve means that it 
has better isolation between the layers than does the CIO4 compound. 

In contrast, the value for [Cu(pz)2(NOs)](PF6) is found 0.13 units higher than the BKT 
curve. Although it has an XY anisotropy smaller by a factor of thirty than those of the CIO4 
and BF4 compounds, it has the highest critical ratio of the three copper/pyrazine compound. 
It is clear that its transition is driven by the 3D interactions. Consequently, its J' / J ratio of 
3 x 10~ 3 , as derived from Eq. (j3J), is a good approximation. [Cu(pz) 2 (N0 3 )](PF 6 ) is therefore 
a remarkably isotropic quasi-Heisenberg copper compound but the least well-isolated of the 
three compounds in this study. 



D. Field-Induced Anisotropy 

In addition to interlayer interactions and intrinsic exchange anisotropy, there is a third 
contribution to the Hamiltonian (Eq. (CQ)) that can induce an ordering transition in a 2D 
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QHAF, the magnetic field. It has been recognized for over a decade 46J] that an external field 
is equivalent to an easy-plane anisotropy in its ability to induce BKT transitions in layered 
Heisenberg antiferromagnets. This work was later extended to show that an external 
field could induce a minimum in the susceptibility of an isotropic (A= 0) 2D S = 1/2 
antiferromagnet. In the case of an XY magnet, the field would enhance the minimum in 
Xz and overcome the intrinsic anisotropy to induce a minimum in \x for sufficiently strong 
fields. 

The effects of field-induced anisotropy on the ordering temperature in a quasi-2D QHAF 
was first reported in the 1995 study of S^CuC^C^ single crystals by NMR and suscepti- 
bility 39). The zero-field Tn for this compound is 256.5 K. When the field was applied in 
the easy plane, Tn increased by 2.3% to 262.5 K in a 4.7 T field, and by 3.4% to 265.3 K in 
an 8.2 T field; when the field was applied normal to the planes, no enhancement in Tn was 
observed. 

More dramatic field-induced effects have been observed in the CIO4, BF4, and com- 
pounds, made possible by their exchange strengths being only ~ 1% as large as that for 
Sr2Cu02Cl2. For external fields exceeding the respective anisotropy fields (Tabled]), minima 
occur for the \ x of each of the compounds. Evidence of this effect is found in Fig. [8] for PF 6 . 
These data were collected in a 1 kOe field, much stronger than the 70 Oe anisotropy field for 
this compound. As revealed by the low-field ac-susceptibility measurements (Fig. [8]) only 
when the dc-field was reduced to less than 70 Oe did the minimum in \ x vanish. In addition, 
the Neel temperatures of all compounds strongly increased in the presence of larger fields. 
For the CIO4 compound, Tn rose from 4.25 K in zero-field to 5.7 K in a 15 T field, an increase 
of 34%. Unexpectedly, the field-induced enhancement of T N was the same independent of 
the orientation of the field. A detailed study of these effects and a discussion of the magnetic 
phase diagram are in progress and will be reported in a subsequent publication. 



V. SUMMARY 



The magnetic properties of three molecular-based quasi-2d S=l/2 Heisenberg antifer- 
romagnets ((Cu(pz) 2 (0104)2, Cu(pz) 2 (BF 4 )2, and (Cu(pz) 2 (C104)2 have been investigated. 
The ordering temperatures have been determined both by /z + SR and susceptibility stud- 
ies. Values of the intralayer and interlayer exchange strengths (J, J') and the exchange- 
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anisotropy parameters A have been determined. Cu(pz) 2 (0104)2 and Cu(pz) 2 (BF 4 )2 have 
similar exchange strengths (17.5 and 15.3 K, respectively), similar critical ratios T-^/J (0.243 
and 0.248, respectively), and similar ratios of anisotropy fields to saturation fields (5.3xl0 -3 
and 5.8xl0~ 3 , respectively). [Cu(pz) 2 (N0 3 )](PF 6 ) has a a weaker exchange strength, a much 
smaller anisotropy field, but the highest critical temperature; it is the least 2D but the most 
Heisenberg-like of the three compounds. The experimental data are all explained by the ex- 
istence of small XY contributions to the spin Hamiltonian, as well as weak coupling between 
the magnetic layers. The influence of field-induced anisotropy on the magnetic behavior has 
been observed. 
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